Constraining Isocurvature Fluctuations with the Planck Surveyor 
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We consider the detection possibilities of isocurvature fluctuations in the future CMB satellite 
experiments MAP and Planck for different cosmological reference models. We present a simultaneous 
10 parameter fit (8 for the case of open model) to determine the correlations between the cosmological 
parameters, including isocurvature cold dark matter contribution to the anisotropy. Assuming that 
polarization information can be fully exploited, we find that an isocurvature perturbation can be 
detected by the Planck Surveyor if the ratio of the initial isocurvature and adiabatic perturbation 
amplitude is larger than 0.07. In the absence of polarization data, the signal from isocurvature 
perturbations can be confused with tensor perturbations or early reionization effects, and the limit 
is larger by almost an order of magnitude. 
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I. INTRODUCTION 

A second generation of satellite experiments, MAP 
and Planck |2j , will soon provide highly detailed temper- 
ature maps of the cosmic microwave background (CMB). 
From the measured spectrum of temperature fluctuations 
one may infer the properties of the cosmic fluid, such as 
its density or the rate by which it comoves with the ex- 
pansion of the Universe, and thus determine most cosmo- 
logical parameters to an accuracy of a few percent ||-|| . 
In addition, the spectrum of temperature fluctuations al- 
lows us to study and constrain physics at very small dis- 
tance scales, corresponding to energies much higher than 
can be achieved in any particle accelerator in the fore- 
seeable future. This is so because the features on the 
temperature map depend on the nature of the primor- 
dial density fluctuation. 

Cosmic inflation is currently the most popular expla- 
nation for the origin of the primordial perturbations, and 
it typically predicts a near scale invariant spectrum with 
Gaussian fluctuations. The perturbations are adiabatic 
with the number density proportional to entropy density 
so that 5(n/s) — 0. This is so because the quantum fluc- 
tuations of the field responsible for inflation, called the 
inflaton, give directly rise to perturbations in the energy 
density of the inflaton field. 

However, the inflaton may not be the only field which 
is subject to quantum fluctuations during inflation. In 
fact, any effectively massless scalar field will fluctuate by 
virtue of the nearly constant energy density of the in- 
flation era with a root mean square amplitude H/(2w), 
where H is the Hubble parameter during inflation. Such 
fluctuations will contribute to the perturbations in the 
microwave background. They may be adiabatic, in which 
case it is difficult to entangle their effect from the infla- 
ton fluctuations. The other possibility is isocurvature 
fluctuations, which are fluctuations in the number rather 



than the energy density so that S(n/s) ^ 0. Examples of 
particle physics models giving rise to isocurvature fluctu- 
ations include axions |||7|], inflation with more than one 
inflaton field and supersymmetric theories with flat 
directions in the potential ||. 

Isocurvature perturbations p~0| , [TT[ ] in a given particle 
species have 8p = with the overdensities balanced by 
perturbations in other particle species, such as radiation. 
At the last scattering surface (LSS) the compensation for 
the isocurvature perturbations can be maintained only 
for scales larger than the horizon, effectively generating 
extra power to photon perturbations at small angular 
multipoles /. As a consequence, the spectrum of isocur- 
vature perturbations differs a great deal from adiabatic 
perturbations, and a purely isocurvature cold dark mat- 
ter (CDM) perturbation spectrum is in fact already ruled 
out on the basis of COBE normalization and erg, 
the amplitude of the rms mass fluctuations in an 8/i _1 
Mpc -1 sphere. (There are suggestions that a decaying 
CDM model could sustain purely isocurvature perturba- 
tions |l3f| ; however, such models are not motivated from 
particle physics point of view). A small isocurvature con- 
tribution might, however, be beneficial in the sense that 
it could help to improve the fit to the power spectrum in 
fio = 1 CDM models with a cosmological constant . 

The simplest possibility is that there is a single CDM 
which has both adiabatic and isocurvature perturbations. 
An example of such a situation is obtained in super- 
symmetric models with the so-called Affleck-Dine (AD) 
fields which also provide a basis for baryogenesis. 
The AD fields are superpositions of squark and slepton 
fields corresponding to the flat directions. During infla- 
tion they will fluctuate along the flat directions forming 
a condensate. They are complex fields and, in the cur- 
rently favoured D-term inflation models jl6| , is effectively 
massless during inflation. The condensate does not cor- 
respond to the state of lowest energy but fragments into 
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non-topological solitons ]17[ , which carry baryonic (and 
sometimes leptonic) number Jlj| and are called B-balls 
(L-balls). The fluctuations in the AD field are both adi- 
abatic and isocurvature |^], and will be inherited by the 
B-balls, which are either stable or decay later into lightest 
supersymmctric particles, which form the CDM. In the 
latter case there will be a lower bound on the amplitude 
of the isocurvature perturbation. 

Adiabatic perturbations |l^-^l| are characterized by 
the gauge-invariant quantity £, which equals Sp/(p + p) 
when the perturbation exits or enters the horizon. Out- 
side the horizon, £ docs not change with time. Isocur- 
vature perturbations jlO|,^0],[n| are characterized by the 
gauge-invariant entropy perturbation S = <5(n c /s 7 ) — 
5 C — (3/4)(5 7 , where n c is the number density of CDM, 
Sj is the entropy density associated with photons, and 
5 C and 5 1 are the relative overdensities in the CDM and 
photon energy densities. During radiation domination 
S ~ 8 C = 5n c /n c . Outside the horizon, S does not 
change with time. "Initial" scale-free (n = 1) adiabatic 
and isocurvature perturbation spectra are given by 
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Following these remarks, we assume in this paper that 
there is just a single species of CDM, and that it has 
both adiabatic and isocurvature fluctuations. As is nat- 
ural when both fluctuations have their origin in infla- 
tion, we also assume that their power spectra are similar. 
The goal of the paper is then to study how well parti- 
cle physics models with isocurvature fluctuations can be 
tested by the future CMB satellite experiments, and at 
what level isocurvature perturbation can be detected. In 
Sec. II we define the Fisher information matrix needed 
for the cosmological parameter error estimates, and de- 
scribe the expected performance of MAP and Planck. In 
Sec. Ill we discuss parametrization of isocurvature fluc- 
tuation and possible cosmological models. In Sec. IV we 
present the results of a numerical study for MAP and 
Planck, while Sec. V contains a discussion of the results. 



II. ACCURACY OF PARAMETER 
DETERMINATION 

For a CMB temperature map with Gaussian fluctua- 
tions, all statistical information is contained in the angu- 
lar power spectrum 



Cti = (\aL\ 2 )- 
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This spectrum contains a wealth of cosmological infor- 
mation. The major cosmological parameters can be de- 
termined from a high-quality measurement of the spec- 
trum, with an unprecedented accuracy. This is because 
the CMB anisotropy reflects conditions at the early uni- 
verse, as well as the geometry out to the present horizon, 



the microwave photons having travelled through essen- 
tially the whole observable universe as well as most of 
its history (2^]. There are however, degeneracies with re- 
spect to some cosmological parameters; different effects 
causing similar features in the power spectrum. Also the 
information from the large angular scales is severely lim- 
ited by cosmic variance. 

Important additional information can be obtained by 
the measurement of the polarization of the CMB [^3] . 
Since polarization of the CMB is caused by scattering 
only, polarization carries information about the LSS only 
(including late scattering due to reionization) , eliminat- 
ing confusion with linc-of-sight effects. Polarization is 
especially useful in separating out the effects of tensor 
fluctuations, isocurvature fluctuations, and reionization, 
all of which produce excess power at large scales. 

Unfortunately the polarization is expected to be of 
rather small amplitude, at least an order of magnitude 
below the anisotropy signal, making it difficult to detect 
p4j . The contamination with the foreground may also 
be a worse problem than for the temperature anisotropy, 
since there are polarized foreground sources, e.g., galac- 
tic synchrotron emission. The MAP satellite is expected 
to provide a definite detection of CMB polarization. The 
higher sensitivity of the Planck Surveyor will then help in 
getting a reasonable power spectrum for the polarization. 

The polarization pattern on the sky can be separated 
into two components, the E-mode and the B-mode [^3| . 
Thus an experiment which also measures the polariza- 
tion, will produce three maps: temperature, E-mode po- 
larization, and B-mode polarization. The B-mode polar- 
ization of the CMB is uncorrelated with the temperature 
or the E-mode polarization. Thus the statistical infor- 
mation is carried by three angular power spectra 
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and one correlation function 



Cei = 
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giving the correlation between the temperature and the 
E-mode polarization p3"|-p7f. The different spectra are 
demonstrated in Fig. ^ 

These four spectra can be calculated for different val- 
ues of the cosmological parameters Sj and compared to 
data from CMB experiments. This way the true values 
of the cosmological parameters can be determined. 

Assuming the actual values of the parameters are Si, 
the information obtained on these parameters by ob- 
servation is described by the Fisher information matrix 

liJUIH 



dCxi 



\ Sz Cov~\Cxi,C Yl ) 



l X.Y 



dCyi 
dsj 
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where C X i = Cti, Cei, Cbi, or C, 
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FIG. 1. The angular power spectra for SCDM (see Table 
II) adiabatic (a) and isocurvature (b) perturbations. These 
scalar perturbations produce E-mode polarization only, thus 
no B-mode spectra are shown. The dashed lines indicate neg- 
ative values (for Cci)- 

The covariance matrices, whose nonzero components 
are 



Cov(C T iC'ti) = m \C T i + 
Cov(CeiC E i) = Ni [c E l + 
Cov{C Bl C B i)=N l [c m + 
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depend on the sensitivity and angular resolution of the 
instrument with 
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where a c is the pixel noise (sensitivity) and 6*f w hm iC is the 
beamwidth (full width at half maximum) of instrument 
channel c. 

The prefactor is given by 



N, = 



(2Z + 1)/ Sky ' 



(9) 



where / s ky takes into account the fact that the whole sky 
cannot be mapped because of foreground contamination. 
We adopt / s k y = 0.65, i.e., assume that foreground can 
be completely removed from 65% of the sky and the rest 
of the sky is unusable. 

The Planck Surveyor carries two instruments, HFI 
and LFI. For MAP and LFI we use the three highest 
frequency channels, for HFI the three lowest frequency 
channels for our estimates (see Table I). 



instrument frequency 9fwhr. 



MAP 



LFI 



HFI 



40 GHz 0.47 

60 GHz 0.35 

90 GHz 0.21 
44 GHz 
70 GHz 
100 GHz 
100 GHz 
143 GHz 
217 GHz 



cxp 

V2 x 35/iK 
V2 x 35ajK 
V2 x 35a*K 
s/2 x 2.4 x 10~ 6 
3.6 x 10~ 6 y/2 x 3.6 x 10" 6 
4.3 x 10~ 6 V2 x 4.3 x 10" 6 



35 mK 
35 

35 iiK 
23' 2.4 x 10" 6 
14' 
10' 

10.7' 1.7 x 10~ 6 
8.0' 2.0 x 10~ 6 
5.5' 4.3 x 10" 6 



3.7 x 10' 
8.9 x 10" 



TABLE I. The beamwidths and sensitivities of the satellite 
experiments |ll 
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The To error on parameter Sj is 



As, 



(10) 



when all parameters are estimated from the same data. 
(If all the other parameters were somehow known a priori, 
the error would be Asi — (Fjj) _ 2). This is a lower limit 
for the error; but when the error is small it is usually a 
good approximation for the error plJH]- 



III. THE PARAMETERS AND MODELS 

We consider a simultaneous fit of 10 parameters, C2, 
Q = ^A + ^m, ^A-^m, Q b h 2 , J] H DM, #0 , t, ns, r, and 
a. Here r is the optical depth to LSS and r = Cj/Cf 
is the tensor/scalar ratio. We assume three matter com- 
ponents, fl m = flcDM + ^hdm + ^b- Note that we have 
defined f2o to include the vacuum energy; it is the param- 
eter which determines the curvature of space. We assume 
T = 2.728 K, F He = 0.24, and N v = 3. For hot dark 
matter we assume all 3 neutrinos have equal masses. 

The isocurvature contribution to anisotropy is given 
by the parameter a, which gives the fraction of due 
to isocurvature perturbations. Thus 



C£°/(C™° + Cf) 



(11) 



Our a is the same as 1 — a of Stompor et al Ej). The 
isocurvature fluctuations are assumed to be primordial 
density fluctuations in CDM, initially balanced by radi- 
ation+baryons+neutrinos to make the total energy den- 
sity homogeneous. The isocurvature spectral index is as- 
sumed to be the same as for adiabatic fluctuations; in 
principle, they could be different, but if both fluctuations 
originate from inflation, and in particular from massless 
field fluctuations during inflation, then their power spec- 
tra are guaranteed to be equal. 

The large-scale anisotropy due to isocurvature perur- 
bations is |2lj 



ST 
~T 



where 



3 3' 



(12) 



(13) 



is the Newtonian gravitational potential |29|. 
whereas for adiabatic perturbations we have just 





fio 




fib 


^HDM 


Ho 


r 


ns 


r 


a 


SCDM 


1.0 





0.05 





50 


0.05 


1 








ACDM 


1.0 


0.7 


0.06 





65 


0.1 


1 








OCDM 


0.4 





0.06 





65 


0.05 


1 









TABLE II. The cosmological models considered. 



ST 1 
— = -<J> 

T 3 ' 



(14) 



where the potential is related to the gauge-invariant £ by 



* = " 2 C 



(15) 



Thus the "observational" isocurvature parameter a is 
related to the amplitudes of the initial isocurvature and 
adiabatic spectra by 



a _ (-6T/T) is 



1 



a 



(-ST/T) 



ad 
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In a given particle physics model A and B are known 
in terms of the model parameters, and Eq. (17) tells how 
these parameters are related to the observed isocurvature 
parameter a. 

We have chosen to consider three different reference 
cosmological models (see Table II), which are the same 
as discussed by Wang et al. p0| : the standard CDM 
model with £Iq = 1, a CDM model with a cosmological 
constant fl\ = 0.7 and f^o = 1, and an open CDM model 
with fl a — 0.4. 

To find the derivatives dC X i/ds t we used CMBFAST 
pl| to run variations around these three reference mod- 
els, with steps Asi or Asi/si between 0.01 and 0.05. 

We COBE-normalized |32] C2 for the three reference 
models. Note that as C2 is one of the parameters to be 
determined from the data, we do not COBE-renormalize 
it for the variations around these reference models. 

Wang et al. (3^] considered 7-parameter fits. Our ad- 
ditional parameters are Uhdm, and a. For OCDM 
we considered 8-parameter fits only, keeping r = and 
Qhdm = 0, since CMBFAST does not allow nonzero val- 
ues for these parameters in open models. 



IV. RESULTS 

We are interested in how accurately the isocurvature 
contribution can be determined and what is the small- 
est contribution that can still be detected. In the fol- 
lowing tables (III, IV, and V) we give the l-o error on 
the isocurvature parameter a for the three experiments, 
MAP, Planck LFI, and Planck HFI, and an ideal exper- 
iment with zero noise and infinitely fine resolution. Full 
sky coverage would reduce these errors by a factor of 

We see that the addition of the two parameters, I^hdm 
and r, increased the error in a determination significantly 
for the case without polarization data, but had only a 
small effect for the case with polarization. This is be- 
cause the temperature spectra from tensor and isocur- 
vature fluctuations resemble each other, resulting in a 
degeneracy between the parameters r and a. (Adding 
Hhdm as a free parameter changes Aa at most by 10%.) 
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Because only tensor fluctuations have B-mode polariza- 
tion, polarization measurements break this degeneracy. 
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FIG. 2. The error ellipses in the (r,a)-plane for (a) SCDM 
and (b) ACDM. The dashed lines are for the case without 
polarization measurement, the solid lines with polarization. 
The four ellipses from the outside inward are for MAP, Planck 
LFI, Planck HFI, and the ideal experiment. 





MAP LFI 


HFI 


ideal 


Without polarization 


SCDM 


0.37 0.26 


0.23 


0.14 


ACDM 


0.38 0.23 


0.20 


0.17 


With polarization 


SCDM 


0.23 0.058 


0.038 


0.015 


ACDM 


0.21 0.067 


0.045 


0.019 



TABLE III. Aa from 10-parameter fits. 



Note that the tensor and isocurvature parameters, r 
and a are non-negative by definition. It is quite likely 
that their true values are smaller than the errors in their 
determination, so that we will only obtain an upper limit. 
For the case without polarization information, there will 
be a degeneracy between these two parameters, such that 
their sum can be obtained with a higher accuracy than 
the parameters themselves (see Fig. |2j). Since the upper 
limit to the sum is also an upper limit to the individual 
parameters, we can then obtain a tighter upper limit to 
a than the Aa given in Table III. This upper limit is 
essentially the Aa given in Table IV, where we have made 
the restriction r = 0. 

For MAP, there are also other significant degeneracies 
among the 10 parameters considered here. For Planck 
without polarization, the only other important degen- 
eracy is with the optical depth r due to reionization. 
Polarization breaks this degeneracy, since early reioniza- 
tion produces a polarization signal at large angular scales 
which cannot be mimicked by other parameter combi- 
nations p^ ]. For Planck with polarization there are no 
significant degeneracies between a and the other param- 
eters (except, of course, a degeneracy between C2 and a), 
and so the errors would not become much smaller even if 
all the other parameters were assumed known (see Table 
V). 



MAP LFI HFI ideal 

Without polarization 
SCDM 035 ail 0T2 0.10 

ACDM 0.22 0.16 0.15 0.15 

OCDM 0.27 0.13 0.064 0.053 



SCDM 0.21 0.056 0.037 0.015 

ACDM 0.16 0.062 0.044 0.018 

OCDM 0.14 0.046 0.033 0.012 



TABLE IV. 


Aa from 8-parameter fits (Rhdm = 0, r 


= 0). 




MAP LFI 


HFI 


ideal 


Without polarization 


SCDM 


0.053 0.049 


0.047 


0.046 


ACDM 


0.055 0.052 


0.051 


0.050 


OCDM 


0.041 0.038 


0.036 


0.035 


With polarization 


SCDM 


0.053 0.039 


0.030 


0.013 


ACDM 


0.054 0.044 


0.036 


0.016 


OCDM 


0.041 0.030 


0.026 


0.010 


TABLE V. 


Aa from 2-parameter fits (C2 and a) 


This 



table is shown to illustrate, by comparison to Table IV, the 
effect of degeneracies with other parameters. Because the best 
determinations of many of the other parameters have to come 
from the same CMB data, a cannot in reality be determined 
with this accuracy. 



.n ? 




With polarization 
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V. DISCUSSION 

From Tables III and IV we find that the (1-cr) detec- 
tion limit of isocurvature fluctuations is a ~ 0.04 for the 
Planck Surveyor. This corresponds to a ratio B/A ~ 0.07 
between the amplitudes of the initial isocurvature and 
adiabatic perturbation spectra. This is at the level pre- 
dicted by certain B-ball models ||, which thus could be 
tested by Planck. It is likely that other particle physics 
models can also be constrained in a significant way on 
the basis (or the absence) of isocurvature fluctuation 
in CMB. However, to reach this sensitivity polarization 
measurements are crucial, as is evident from Tables III 
and IV. 

Without polarization information, the contribu- 
tion from isocurvature fluctuations to temperature 
anisotropies could be confused with tensor perturbations 
or early reionization effects. Polarization measurements 
break this degeneracy, since these other two effects both 
have a unique polarization signal. 

Whether such sensitivity can be realized is still un- 
clear; the foreground contamination for the polarization 
spectra may turn out to be a more severe problem than 
in the case of the temperature angular power spectrum. 
However, optimizing the polarization information should 
have a high priority, in particular for Planck; because of 
its higher sensitivity, it is in this respect that Planck is 
truly a superior instrument compared with MAP. CMB 
polarization measures both tensor and isocurvature con- 
tribution, which are of paramount interest to inflation 
model builders. 

Finally, the cosmic variance sets an absolute lower 
theoretical limit, at about a ~ 0.01, to how small 
isocurvature fluctuations can be detected from the CMB 
anisotropy by any experiment. 
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